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Abstract 

This note addresses a three-dimensional model for isothermal stress-induced transformation 
in shape-memory polycrystalline materials. We treat the problem within the framework of the 
energetic formulation of rate-independent processes and investigate existence and continuous 
dependence issues at both the constitutive relation and quasi-static evolution level. More- 
over, we focus on time and space approximation as well as on regularization and parameter 
asymptotics. 
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1 Introduction 

Shape-memory materials are metallic alloys showing some surprising thermo-mechanical behavior: 
severely deformed specimens with residual strain up to 15% regain their original shape after a 
thermal cycle {shape- memory effect). Moreover, the same materials are super-elastic (also called 

*This research was partially supported by the European Union via HPRN-CT-2002-00284 Smart Systems: New 
Materials, Adaptive Systems and their Nonlinearities and by the DFG Research Center 86 Matheon under sub- 
project CIS. 
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pseudo-elastic) , namely, they recover comparably large deformations during mechanical loading- 
unloading cycles at prescribed temperatures (see, among others, [TJ [B] I^TJ [23 SSI HZ] ) ■ These 
features, which are not present (at least to this extent) in materials traditionally used in engineer- 
ing, are at the basis of the innovative and commercially valuable applications of shape-memory 
materials. Namely, shape-memory technologies are nowadays exploited in a variety of different ap- 
plicative contexts ranging from sensors and actuators (even microscopical), to robotics, to clamp- 
ing and fixation devices, to space applications (grippers, positioners), to damping devices (shock 
absorption) [56]. The largest commercial success of shape- memory materials is however related 
to biomedical applications. The combination of good bio-compatibility and interesting material 
properties creates unique materials for medical tools and devices. Nowadays, shape-memory mate- 
rials are successfully used in orthodontics (archwires), orthopedics (bone anchors, intromeduUary 
fixations, bone staples), medical instruments, minimal invasive surgery technology (catheters, en- 
doguidewires, grippers, cutters), drug delivery systems, and both intravascular (cardiovascular 
stenting, bronchial biliary, aortic aneurysm, carotid stenosis) and extravascular scaffolding. In 
particular, shape-memory stents are the key tool in order to implement a variety of quite success- 
ful non- invasive surgical techniques [HI [Ml [55] . 

The present analysis is concerned with the quasi-static evolution of shape-memory materials in 
the small-strain regime. In particular, we shall study a macroscopic phcnomcnological model for 
shape-memory polycrystalline materials undergoing stress-induced transformations that was origi- 
nally proposed by SouzA et AL. [53] and later addressed and extended by AuRiCCHiO & Petrini 
[4l[5], and AuRiCCHiO et AL. [7|. Our aim is to focus on the isothermal situation at suitably high 
temperatures in order to capture the super-elastic material behavior. The understanding and the 
efficient description of the super-elastic regime is clearly of a great applicative interest. In partic- 
ular, most of the biomedical applications enlisted above are based on super-clastic deployment in 
situ and/or super-elastic kink resistance of shape-memory materials. 

Let us briefly recall here the basic features of the proposed model, the interested reader is 
of course referred to the above-mentioned contributions for all the necessary modeling details 
and motivations as well as for some computations and validation. The formal character of this 
introduction is intended to serve for the purpose of a general overview on the model and our 
results. In particular, (most of) the mathematical details are here omitted and will be provided in 
the forthcoming sections. 

Moving into the frame of Generalized Standard Materials (see Maugin ^9\) and within the 
small-strain regime, we additively decompose the linearized deformation e = (ffy ) = (uij + Uj,i)/2, 
(u being the displacement from a fixed reference configuration fl C M.^) into the clastic part Sei 
and the inelastic (or transformation) part z as 



At the microscopic level the super-elastic effect is interpreted as the result of a structural phase tran- 
sition between different configurations of the material lattices, namely the parent phase (austenite 
and twinned martensite) and its shared counterpart termed product phase (detwinned martensite). 
In particular, the internal variable z is assumed to be descriptive of the mechanical (tensorial) 
effect of the detwinning observed in the material. 

Denoting by W{e^ z) the stored energy density of the system, the evolution of the material will 
be described by the following classical relations 



£ = £cl + Z. 



(I.I) 



z 



dW/de, 
dW/dz, 



(1.2) 
(1.3) 
(1.4) 



Here, ^ denotes the thermodynamic force associated with z and (|I.4p is the flow rule for z where 
D* stands for the Legendre conjugate of the dissipation density D (see below). 
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The material constitutive relations p.2p - p.4[) may be conveniently rewritten in the following 
equivalent subdifFerential formulation 



\ fd,Wie,z)\ fa 
dD{z))^\d,W{e,z)) \0 



(1.5) 



where D stands for the dissipation density and the symbol d denotes subdifFerentials in the sense 
of Convex Analysis (see below). 

The evolution problem (jl.Sp may be set within the frame of energetic formulations of rate- 
independent processes recently proposed by Mielke et AL. [28l [42l [44] . The notion of energetic 
solution (discussed in some detail in the forthcoming Section [2]) is based on equivalently recasting 
the subdifferential problem p.5p as the coupling of a global stability condition and an energy 
conservation relation. In particular, the subdifferential relation (jl.Sp is rewritten as 

(stability) {e{t),z{t)) eAi-gMin(w{e,z)~a{t) -.e + Diz-zit))) (1.6) 

ie.z) ^ ' 

(energy equality) W{e{t), z{t)) - a{t) : e{t) + DissD(z, [0, t]) 

= W{eo,zo)-a{0):eo- f a{s) : e{s) ds, (1.7) 



for all t > 0. Here, we assume to be given some suitable initial data (eoj^o) and the stress 
1 1-^ a{t) and denote the total dissipation of the system on [0,t] as 



Dissr,(z, [0, t\) := sup <^ V D{z{U) - z(ij_i)) : {0 = to < ii < • • • < ^at-i < tN = t} 




where the supremum is taken with respect to all finite partitions of [0,t\. Energetic formulations 
were originally developed for shape- memory alloys in Mielke & Theil and Mielke et al. 
j42[ l43l [44] , and have shown to be extremely well-suited for a variety of different rate-independent 
situations. In particular, they have been successfully considered in connection with elasto-plasticity 
[ni[Tl[31[31[3Sl[SS], damage [3S], brittle fractures [H], delamination gS], ferro-electricity [i5] . 
shape- memory alloys |38[|42[[^ . and vortex pinning in superconductors [52| . The reader is referred 
to Mielke (23 for a comprehensive survey of the mathematical theory. 

Let us now introduce the precise form of W we will deal with. Namely, we choose 

Wie,z) - ic(e-z) : {e - z) + ci\z\ + C2\z\^ + I{z) + ^\Vz\\ (1.8) 

Here, C is the elasticity tensor and the positive parameters ci and C2 are given. Indeed, in 
[53] the constant ci is assumed to depend explicitly on the temperature of the specimen while 
here temperature effects are neglected. On the other hand, C2 measures the occurrence of some 
hardening phenomenon with respect to the internal variable z. The function / is the indicator 
of a fixed closed ball of radius C3 > 0. In particular, C3 represents the maximum modulus of 
transformation strain that can be obtained by alignment (detwinning) of the martensitic variants. 
Finally, the positive coefficient v is expected to measure some nonlocal interaction effect for the 
internal variable z and Vz stands for the usual gradient with respect to to spatial variables. 
Indeed, gradients of inelastic strains have already been considered in the frame of shape-memory 
materials by Fremond [20] and the reader is referred also to Arndt et al. |2], Fried & 
GuRTiN [22], Kruzik et AL. [2S], MiELKE & RoUBi'cEK [3H] , RouBiCEK [SHI [S] for examples 
and discussions on nonlocal energy contributions of z. 

The proposed model is capable of describing the main features of the super-elastic evolution of 
shape-memory materials. In particular, the internal variable tensorial character of the model allows 
for taking into account the so-called single-variant martensite reorientation phenomenon. Namely, 
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also in the case the material is fully transformed into product phase (i.e. \z\ = C3), inelastic strain 
changes can still be experienced due to variant reorientation (i 7^ 0). This fact is experimentally 
observed and turns out to be crucial with respect to applications. Moreover, whenever not restricted 
to the isothermal situation, the model turns out the be thermodynamically consistent in the sense 
that the Second Law of Thermodynamics is satisfied in the form of the Clausius-Duhem inequality. 

As for the full quasi-static evolution of the material we shall couple the constitutive relation 
(jl.Sp with the equilibrium equation 

divcr + / = in 17, (1.9) 

where / is a given body force, suitably complemented with some prescribed boundary displacement 
and boundary traction in distinguished parts of the boundary of 17. 

The first issue of this paper is that of adapting the above referred abstract theory for energetic 
formulations to the quasi-static evolution problem and obtain that fTheorem lG.ip 

(existence) the quasi-static problem admits at least one energetic solution t ^ {u{t),z{t)). 

We shall be concerned with some specific rcgularization of the original quasi-static model. 
Namely, some smooth variant of the potential W above turns out to be better suited for the sake 
of numerical considerations. In particular, we will consider a regularized version of the model by 
posing 

W^p,,(£, z) = ic(£ -z):{e-z) + Fp{z) + ^jVzj^, (1.10) 

where > and Fp is some regularization of Fq : z ^ ci\z\+C2\z\'^+I{z) obtained by penalization 
and smoothing and depending on the regularization parameter p > 0. This regularization is exactly 
the starting point of AuRiCCHiO & Petrini [HE], and has been exploited in AuRiCCHiO et AL. 
[7] as well (in all these papers 1/ = though). 

A second focus of the present contribution is on unique solvability of the regularized model. In 
particular, we check that 

(uniqueness for p > 0) for p > 0, the quasi-static problem has a unique solution. 

This uniqueness result was proved in an abstract frame by Mielke & Theil [1^ H5] and is here 
reconsidered in the specific situation of the regularized version of the quasi-static problem. 

A quite natural approach to rate-independent evolution problems relies on implicit time- 
discretization. This perspective is here investigated and complemented with some space approx- 
imation technique. In particular, the main novelty of this paper is the convergence analysis for 
the discretized-regularized model. Namely, we consider the (possibly joint) limits with respect 
to the time-steps r of time partitions (here considered to be constant for simplicity), the space 
mesh size h (conforming finite elements are exploited), and the regularization parameter p. In 
particular, denoting by (m, z)p^T,h the unique solution to the space-time discrete problem with 
the parameter-choice p > (time-interpolant, piecewise constant on the time-partition) and by 
(m, z)p the time-continuous solution to the problem for p > 0, we prove the following (Theorem 

(convergence for p > 0) for p > 0, {u,z)p^r,h converges to {u,z)p as (t, /i) — ^ (0,0), 
(full convergence) up to a subsequence, {u, z)p^r,h ^ {u, z)o as (p, r, ft,) (0, 0, 0). 

Of course the topologies under which the latter convergences hold true will be specified in the 
forthcoming sections. 

Indeed much more is true and we are in the position of giving a full picture of convergences 
for the model subsequently. Moving from Section [5] where the mathematical formulation of the 
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problem is presented, we shall organize our results by successively increasing complexity. Section 
[3] addresses the analysis of the constitutive relation problem (jl.5|) . namely the zero-dimensional 
problem. In particular, we prove well-posedness and convergence of time-discrete approximations. 
Then, the three-dimensional minimum problem arising from time-discretization is addressed in 
Section [4] where we also investigate well-posedness and convergence of space approximations along 
with suitable error bounds. Some a priori bounds and a preliminary convergence result for the 
incremental solutions to the problem in case the time-partition is fixed are discussed in Section 
[5l Finally, the three-dimensional quasi-static evolution problem is tackled in Section [6] where we 
provide the above mentioned existence, uniqueness, and convergence results for the space-time 
discrete solutions. Finally, Section [7] deals with convergence issues with respect to parameters and 
discretizations in full generality. 

2 Mathematical formulation 

Tensors. We will denote by M.'l^^ the space of symmetric 3x3 tensors endowed with the 
natural scalar product a : b :^ tr(a&) = aijbij (summation convention) and the corresponding 
norm |ap := a : a for all a, b £ Mgy^. The space Mgy^-^ is orthogonally decomposed as ^^sym — 
■'^dev ffiRl2, where ]Rl2 is the subspace spanned by the identity 2-tensor I2 and K'^ov t'^^ 
subspace of deviatoric symmetric 3x3 tensors. In particular, for all a e Kgymi have that 
a = fldev + tr(a)l2/3. For aU u £ H^^^{R^;R^) we let e{u) e Lf^^{R^;Rl^l) denote the standard 
symmetric gradient. 

Reference configuration. We shall assume 51 to be a non-empty, bounded, and connected 
open set in with a Lipschitz continuous boundary. The space dimension 3 plays essentially 
no role throughout the analysis and we would be in the position of reformulating our results in 
R'' with no particular intricacy. We assume that the boundary dfl is partitioned in two disjoint 
open sets Ft,- and Foir with 9Ftr = 9FDir (m We ask Foir to be such that there exists a 

positive constant cq depending on FDir and f2 such that the Korn inequality 

co||w||ffi(n;R3) < ll"lli2(rD..;K3) + \\^(^)\\l2(^n,Ml^i)^ (2-1) 

holds true for all u g H^{n; R'^). It would indeed suffice to impose Foir to have a positive surface 
measure (see, e.g., HH Thm. 3.1, p. 110]). 

Prescribed boundary displacement. We will prescribe some non-homogeneous Dirichlet 
boundary conditions on FDir- To this end, we will assign u^" G C^([0, T]; i?^/^(FDir, R'^)) or, 
equivalently, vP" G C^{[0,T]; H^{i},R^)) whose trace on Foir is the prescribed boundary value 
for the displacement u. On Ftr some time-dependent traction will be prescribed instead. 

Elastic energy. Let C be the elasticity tensor. The latter is regarded as a symmetric positive 
definite linear map C : Rsy,Jj I^sym- We shall assume that the orthogonal subspaces R'^^.^ and 
RI2 are invariant under C. This amounts to say that indeed 

Ca = CdovOdov + Ktr(a)l2, 

for a given Cdcv : H^dcv ~^ ^dov ^^"^ ^ constant k, and all a G Rgym- The case of isotropic 
materials is given by Cdcv = 2G(l4 — I2 ^ I2/3) and G and k are respectively the shear and the 
bulk moduli. The latter decomposition is not exploited in our analysis but it is clearly suggested 
by the mechanical application. 
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We will make use of the stored elastic energy functional C : L'^{fl] Kgym) ~^ [Oj +00} defined as 

2 Jo 



Inelastic energy. As for the stored inelastic (or transformation) energy we shall prescribe the 
function F : M^^^^ ^ [0, +00] as 

F{a) = ci\a\ + C2\af + lia), 

where / : M^^^^ ^ [0, +00] is the indicator function of the ball {a G R^cv ■ l*^! — "^si ^^'^ ^^"^ 
positive constants ci, C2, and C3 are given. Moreover, the stored inelastic energy functional is 
defined as T : L'^{n]R^^J) [0, +00] as 

T{a) / F{a)dx if F{a) G L^{n) and T{a) = +00 otherwise. 
Jn 

The well-posedness and time discretization issues discussed here do not rely on the particular 
form of F and could be adapted to any uniformly convex, proper, and lower semicontinuous 
function. We however prefer to stick to the actual modeling choice for the sake of clarity. In the 
forthcoming of the paper we will address some suitable regularization of F. Indeed, we introduce 
an approximation parameter p > and some functions 

FpGC^^^Rf^^) with V^Fp bounded, V^Fp > C2I4, and Fp(0) = 0, (2.2) 

and define Fo := F. An example in the direction of (|2.2[1 is 

Fp{a) ci(v/p^ + |a|2 - p) + C2\af + ip{\a\)/p 

for ipeC^'\R), ip'eL°°{R), ip">0, </j(r) = iff r < C3. (2.3) 

Exactly as above, for all p > we let the regularized stored inelastic energy functional Tp : 
i2(fj;M3^^3^ [0,+cx)) be defined as 

Tp{a) / Fp{a)dx, 
Jn 

and J-Q := J-. Finally, we shall be considering also some space-regularized situation. To this end, 
let p,iy>0 and define Tp^^ : L'^{n]M.\''^^) [0, +00] as 

TpAa) := [Fp{a) + ^|Vap) dx, 

where {'Va)ijk = daij/dxk is the usual gradient in the distributional sense and | • | denotes here 
the Euclidean norm. 



Stored energy. Following the above introductory discussion, we define the stored (Helmholtz 
free) energy functional for p, v >0 as 

>Vp,^(u,z) ■.= C{e{u)- z)+Tp,^{z). 



Load and traction. We assume to be given the body force / G iy^'^(0, T; M^)) and 

a surface traction g G iy^'^(0, T; L^(rtr; K'^)). In particular, one can define the total load t G 
W'^-^{Q,T; (i7i(f7;M3))') (the prime denotes here the dual) as 

{i{t),u) -.^ [ f-udx+ [ g-udH^ Vm G i/^(f7; M^), t G [0, T], 

where Ti.^ is the 2-dimensional Hausdorff measure and (•, •) denotes the duality pairing between 
(ffi(rj;R3))' and H^{n;R^). 
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State space. We set our problem by letting 

Here = for = and = 1 otherwise. For all u e H\n; R^), let us define ^{11) C y 



as 



yiu) := {{u,z) : u = u on Ton}, 

Then, for all t e [0, T], we shall define the phase space of the process as y'^{u^"{t)). For the sake 



of later purposes (see also (jl-Sp ) let us denote by Wp : Mgy^ x M^^^^ [0, +oo) the function 



Wp{e, z) := ic(£ - z) : (e - z) + Fp{z), 

and remark that, owing to (|2.2p . the second derivative V'^Wp is bounded and Lipschitz continuous. 
Moreover, we denote by Au ■ y ^ [0, oo) the quadratic form 

|2j^, , / |V7-,|2, 



:= C(e(w) - z) + C2 / |z|^da; + - / \Vz\^dx V(w, z) e 
and by a > the corresponding uniform ellipticity constant (depending on C, C2, and v). 



Dissipation potential. The quasi-static evolution of the material is described by means of an 
appropriate dissipation mechanism, see (jl.Sp . To this aim, we choose the dissipation (pseudo)- 
potential D : K'^^^ [0, +oo) to be lower semi-continuous, positively 1— homogeneous, and to 
fulfill the triangle inequality 

D{a) < D{b) + D{c) whenever a = 6 + c. (2.4) 

Moreover, we ask for some constant cq > Q such that 

CD\a\<D{a) VaGM^dcv- 

Under the current assumptions on Z?, the latter non-degeneracy condition is indeed equivalent to 
the fact that the set {a : D{a) < 1} is bounded or that D does not vanish except in 0. Let us 
stress that D turns out to be convex (see (|2.4p ) and that there exists a second constant Co > 
such that 



D{a)<CD\a\ Vae™'''^^ 



dcv ■ 

3x3\ 



We define the corresponding dissipation fimctional V : {Q;Mr^^) [0, -l-oo) as 

V{a) = / D{a)dx. 



One shall stress that indeed, since 2? is obviously positively 1-homogeneous, a rate-independent 
evolution follows. Moreover, we recall here that, for all z : [0,T] — > K^g^, we let 



Dissr,(z, [s,t]) := sup |^ i:»(z(t,) - z(t,_i)) : {s = io < < • • • < ^w-i < = i}| 



(2.5) 



the supremum being chosen on the set of all finite partitions of [s, t] C [0, T]. Finally the analogous 



notion Dissx)(z, [s,t]) will be used for functions which take values in L^(r2; M'^g^). 
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State space approximation. Henceforth wc will be interested in some space approximation 
procedure. Indeed, we assume to be given a suitable sequence of approximating closed subspaces 
3^^ := Uh X C depending on some parameter ft, > which is intended to go to zero in 
the limit. We shall collect and comment here the abstract assumptions which will be exploited in 
the following. Of course the main application we have in mind are conforming finite elements on a 
shape regular and quasi-optimal mesh [8] with size h on the polyhedral domain f2. We will firstly 
ask to be non-decreasing and such that U/i>oy/^ is dense in . Moreover, we restrict from 
the very beginning to the special case when = y\ C 3^^. 

Now let p'^'.y'^ ^ y^ the Galerkin projector corresponding to the scalar product induced by 
the quadratic form Au- In particular, by introducing the bilinear form '-y^ 'X-y^ defined 

by 

B^{{ui,zi), {u2, Z2)) := ^ / C(e(Mi) - zi) : {e{u2) - 22) + C2 / 21^2 + ^ / Vzi • Vz2 

for (ui,zi), (u2,22) £ , we have that, for all {u,z) G , the projection p^{u,z) may be 
uniquely determined by 

B„{[u,z) -p1^{u,z),{uh,Zh)) ^0 'i{uh,Zh) eyi- (2.6) 

Namely, one has that 

A,,{pl{u,z)) = B,{pl{u,z),pl{u,z))<A,{u,z) VKz)er. (2.7) 

Let us explicitly observe that p^ is pointwise converging in to the identity as /i — > 0. 

Next, let us introduce a pair of operators qt '■ U ^ Uh and : Z'^ Z'j!^ and ask them to be 
pointwise converging to the identity as /i — * 0. More specifically, we will ask for 

h^Q,v^Q => rl{z)^z Mz^Z". 

Moreover, we require that 

z G Z" and \z\ < C3 a.e. in f2 =^ k/^l^)! ^ ^3 a.e. in 57, (2.8) 

and that r° : ^ maps bounded sets into bounded sets. As for an example of operator 
fulfilling the assumptions is the component-wise Clement interpolant from L^(r2;R^^^) to the 
space of piecewise linear functions [9]. In this case, relation (|2.8p follows from Jensen's inequality. 



3 Analysis of the constitutive relation 

Let us start our analysis by focusing on the constitutive material relation. Namely, we neglect 
for the moment the coupling of the material model with the equilibrium problem p.9p . Assuming 
to be given a tension history, we solve for the elastic and the inelastic strain starting from a 
given state. The understanding of this simplified (reduced) problem will be crucial. First of all, 
a detailed study of the constitutive relation is surely an important step in the direction of the 
investigation of the full quasi-static evolution problem. This in especially true with respect to 
numerics. Indeed, the efficient solution of the constitutive relation is the key ingredient for a full 
discretization procedure. Secondly, the full equilibrium system might reduce to a zero-dimensional 
problem under specific yet common geometric restrictions or symmetries. Finally, we aim to give 
in this somehow (notationally) simplified situation the main points of our analysis. 

Assuming to be given a : [0, T] Kgymj we shall determine e : [0, T] Kgym and z : [0, T] 
-'^dev starting from (eo,zo) and fuffiUing (jl.Sp . Of course, since the transformation strain z is 
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assumed to be deviatoric and the elasticity tensor C decomposes as above, the problem could be 
easily reformulated in the deviatoric subspace R'^p^ only. We however prefer not to exploit this 
simplification for the sake of consistency with the forthcoming analysis. 

Let p > be fixed throughout this section. We shall be concerned with the energy function 
Wp(e, z) - a{t) : e which is defined for all [t, e, z) £ [0, T] x Rl^^ x R^^^^. Moreover, let us define 
the set of stable states at time t £ [0, T] as 

Sit) := {(£, Z) G R3^3 ^ ^3X3 ^^^^ ^1^^^^ g ^3X3 ^ j^3^x3^ 

Wp{e,z)-a{t) : e < Wp{e,z)-a{t) : £ + (3.1) 

and S ■.= Ute[o,T]{t,S{t)). 

As for an energetic solution of (|1.5p we mean a pair (e, z) : [0, T] Rgym ^ -"^dev such that 
the function t cr(t) : e{t) is integrable and, for all t G [0,T], 

is{t),z{t))eS{t), (3.2) 
Wp{e{t), z{t)) - a{t) : e{t) + Dissi)(z, [0, t]) 

= Wp{eQ,zo) - a{0) : eo - a{s):e{s)ds. (3.3) 

Let us now comment on the equivalence between ()1.5p and the energetic formulation (|3.2p - 
()3.3p . To this end we will focus for simplicity on the smooth case p > 0. Indeed, the argument 
for the situation p = is just slightly less straightforward from a notational viewpoint. Using the 
definition of the subdiffercntial dD{z), relation p.5p turns out to be equivalent to 

(deWpie, z) ~a) : (v ~e) + d,Wp{e, z) : (w - i) + D{w) - D{z) > 

\/iv,w) G R3^1 X Rf^^, a.e. in (0,T). (3.4) 

Now, by respectively choosing (v,w) = {kv,kW) and letting k — > +oo or {v,w) = (0,0) in the 
latter relation we easily get that 

(deWpie, z)-a):v + d,Wp{e, z):w + D{w) > 

y{v,w) G R%1 X Rf^'^, a.e. in (G,T), (3.5) 
{dMp{s,z) ~ a) : s + d,Wp{s,z) : z + D{z) <0 a.e. in(0,r). (3.6) 

Of course p.4p and (|3.5p - (|3.6p are equivalent. Now, since Wp is strictly convex, we have that 
(e(t), z{t)) is the almost everywhere unique minimizer of 

(e,z) 1-^ Wp{e,z) - a ■.e + D{z-z[t)). 

In particular, by assuming e, z, and cr to be absolutely continuous (see below), we readily check 
that dSSl) holds. Moreover imply that 

{d^Wp{e,z)-a):e + d^Wp{£,z):z + D{z) = Q a.e. in (0,T), 

which can be rewritten as 

^(Wp(£,z) -CT : e) = -(T : £-D(i) a.e. in (0,r). 

Hence, by integrating the latter on (0, i) for t G [0, T], we readily deduce (|3.3p . Vice versa, (|3.3p 
allows us to recover p.5p - p.6p at once by differentiating and exploiting (j3.2p . 

The main advantage of the energetic formulation p.2p - p.3p is that it does involve neither 
derivatives of constitutive quantities nor of the solution. It is hence particularly well-suited for the 
aim of proving wcU-posedness results and it simply generalizes to possibly non-convex situations. 

The aim of this section is to exploit here the abstract existence theory for energetic formulations 
developed in [T71 US] and adapt it to the current modeling situation. 
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The incremental problem. In order to find an energetic solution to (|3.2|) - (|3.3p we shall consider 
an implicit time discretization procedure. At first, let us observe that, for all z S I^dov 
t £ [0, T], the function (e, z) i— > Wp{e, z) — a{t) : e + D{z — z) has a unique minimum since it is 
uniformly convex and coercive. Let now the partition P := {0 = to < ti < ■ ■ ■ < t^-i < = T} 
be given with diameter r = maxi=i_...^jv — ^i-i- Moreover, let (eo,^o) € S{0) be a given initial 
datum. One should consider that, for any given zq G K^gv' there exists a unique eg = jCzq, where 
C = id here, with {eq, zq) g <S'(0). Hence, we solve iteratively the minimum problem 

{e,,z,)€ ArgMin {Wp{e, z) - a{ti) : e + D{z - z,.^)) for 2 = 1, . . . , iV. (3.7) 

We shall refer to the latter as the incremental problem associated with p.2p - p.3p . Let us explicitly 
observe that, by the triangle inequality, any solution ) to (|3.7p solves also 

(£„z,)e ArgMin {Wp{e, z) - a{ti) : e + D{z - z,)) for z = 1, . . . , TV. (3.8) 



Error propagation. We shall start by providing a continuous dependence result for the single- 
step minimum problem in (|3.7p . Referring to the forthcoming time-stepping procedure, the follow- 
ing estimate can be seen as some error propagation control. 

Lemma 3.1 (Continuous dependence). Let {a^,z^) G Kgym ^ ■'^dov J ~ 1'^, be given and 



(e^z^) := ArgMin(^_^)gjj3><3^jj3^x3(H/p(e,z) : e + D{z ~ z^)). 



Then 



\S' - + |zl - z2|2 < J_|^l _ a2|2 ^ _ _2 



a 



Proof. Since (£^,z^) is minimal and Wp is uniformly convex of constant a one has that 
a|£i-f2|2 + a|^i_^2|2 < Wp{e^z^)-a' -.e^ + Diz^-z') 

- Wp{e\z^)+a^ -.e^ ~D{z^ ^z'^). 

On the other hand, the minimality of (e^, z^) entails that 

< Wp{e\z^) -a^ -.e^ + D{z^ - z^) - Wp{e^,z^) ct^ : - D{z^ - z^). 
Taking the sum of the latter relations and exploiting the triangle inequality (|2.4p we get that 

a\e^ - e2|2 + a\z^ - zT < i^' - ■ (e' " e') + ^D{zi - za), 
whence the assertion follows. □ 



The evolution problem. We shall now provide the main result of this section which follows by 
passing to the limit in the above described time-discrete approximation. 

Theorem 3.2 (Existence for p > 0). Given a G W^-^{Q,T;W^'^ and (eo, zo) S 5(0) there 
exists an energetic solution (e, z) to (j3.2[l - (j3.3|) such that (e(0),z(0)) = {eo,zo). Moreover (e, z) G 

Proof. Let us choose a sequence of partitions P" := {0 = tg < t" < • • • < t'^„_^ < = T} 
with diameters r" = maxi=i....^jvn (i" — t2_i) going to zero. Owing to the above discussion, 
we uniquely determine a sequence of solutions {(e",z")}^Q to the corresponding incremental 
problems (j3.7p such that (eq , Zq ) = (eg, zg). We shall denote by (e", z") the incremental solution^ 
i.e. the right-continuous piecewise- const ant interpolant of {(e", ^Dlil^o on the partition P", and 
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by r",s" : [0, T] ^ [0, T] the functions r"(0 := for t e and s" (i) := for 

Since {(e^^DliTo solves dSH) with z," replacing Zi, one directly gets that (e^zf) G 5'(ir) 
for all 1 = 1,..., A^". Moreover, from (|3.7p and the minimality of (e", z"), we compute that 

Next, taking the sum of the latter relation for i = 1, . . . , m and m < A^", we get that 
M^p(er„, <J - ^(C) : - Wpie„,z„) + a(0) : Eq 

+ ^ D(zr - zr_i) < - / a : e" ds. (3.10) 

^=l 

Hence, it suffices to apply the discrete Gronwall lemma and exploit the cocrcivity of Wp in order 
to check that 

sup lVp(£"(i), z"(i)) and Diss£)(z", [0, T]) are bounded independently of n. (3.11) 
te[o,T] 

Indeed, the latter bound depends on Wp{eo-,zo) and HcH^yi.i^Q y.jjsxsj only. 

In order to pass to the limit with n we exploit Helly's selection principle and find a (not 
relabeled) subsequence of partitions and a non-decreasing function (f) : [0,T] — > [0,+oo) such that 

z"(t) z(i), Dissi5(z", [0,t]) ^ </-(i) for all t e [0,T], (3.12) 
and Dissi5(z, [s,t]) < 0(t) -0(s) V[s, t] C [0, T]. (3.13) 

Consequently, for all t G [0,r], we readily find the unique limit e{t) = Cz{t) since e"(t) = 
Cznit) Cz{t). 

Next, we check that S is closed. Indeed, let the sequence {tk,ek,Zk) £ S converge to (i,e,z) 
in [0, T] X Mgyi'ji X K^ct- Then, since Wp is lower semicontinuous and a is continuous, for all 

{e,z)eRl^lxRlf^, 

Wp{e, z) - ait) : e < liminf (Wp(efc, Zk) - cr(ifc) : £fc) 

fc— >-t-oo 

< liminf (VFp(e,z) - a{tk) : s + D{z - Zk)) = Wp(e,z) - a{t) -.e + Diz- z). 

k—^+oo 

Namely {t, e, z) £ S. We shall exploit the latter closure property in order to prove that (e(i), z{t)) 
is a stable state. Indeed, recalling that t E [0, T] is fixed, one readily checks that the sequence 
T"(t) converges to t and is such that (e"(T"(t)), z"'(T"'(t))) converges to (e(i),z(<)) by definition. 
Hence, relation (|3.2[) follows since (T"'(i), e"'(T"'(t)), z"'(T"(t))) S 5. In particular, we have proved 
that {e{t),z{t)) solves (see dSll)) 

(e(t),z(0)G ArgMin (Wp(£, z) - cr(t) : £ + L>(z - z(t))) . 

Moreover, by construction, we have (£(0), z(0)) = (£o, zq). 

We are left to prove that indeed (£, z) fulfills the energy identity (|3.3p . Relation (|3.10p can be 
rewritten as 

Wp{e"{t),z"{t)) - (7(t"(<)) : £"(0 +Dissz,(z", [0,r"(0]) 

< Wpieo, zo) - <J{Q) : £o - / a : £" ds. (3.14) 

Jo 
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Hence, passing to the liminf in the latter relation and exploiting once again the lower semicon- 
tinuity of Wp, the integrability of a, the boundedness of e" (see p.ll|) '). and (|3.13p . we readily 
check by Lebesgue dominated convergence that 

Wp{e{t), z{t)) - <j{t) : e{t) + Bissoiz, [0, t]) 

< Wp(£o,zo)-fT(0):eo- / a-.eds. (3.15) 

Jo 

Some more precise convergence for the energy can be deduced. Indeed, from the stability 
condition (e"(t), z" (t)) G S'(s"(t)), the lower semicontinuity of Wp, and the continuity of a one 
checks that 

Wpieit), z{t)) - ait) : e{t) = lim (Wpieit), z{t)) - a(s"(t)) : e(t) + D{zit) - z"(t))) 

>limsup(W^p(£"(i),2"(i)) -f^(s"(i)) : > Wp{e{t), z{t)) ~ a{t) : e{t). (3.16) 

n — S-+00 

In particular, we have proved that Wp(e"(t), z" {t)) converges to Wp{e{t), z{t)). 

Our next step will be that of proving that (e, z) is absolutely continuous. Indeed this follows 
at once from the stability condition (|3.2p . the upper energy estimate p.l5|) . the uniform convexity 
of Wp, and the absolute continuity of a. Let us fix [s,t] C [0,r]. Owing to (e(s),z(s)) £ S{s) 
and the uniform convexity of Wp with constant a one readily gets that 

a\e{t) ~ e{s)\^ + a\z{t) ^ z{s)\^ 

< Wp{e{t), zit)) - ais) : e{t) + D{z{t) ~ z{s)) ~ Wp{e{s), z{s)) + a{s) : e{s) 

< Wpieit), zit)) - ait) : e(t) + Dissi3(z, [s, t\) 
-Wpieis), zis)) + ais) : £(s) - (a(s) - ait)) : sit) 

< - / air) : (e(r) -e(t))dr. 



Hence, by means of Gronwall's lemma, one checks that 

\eit) - eis)\ + \zit) ~ zis)\ < a I \a\, (3.17) 



where the positive constant C4 depends just on a. The absolute continuity of e and z follows. 

We are now in the position of proving the converse inequality to (|3.15p . namely, the lower 
energy estimate. Indeed, for all t S [0,T], 

Wpieit), zit)) - ait) : £(t) + Bissoiz, [0, t]) 

>Wpieo,zo)-aiO):ea-[ a:eds. (3.18) 

Indeed, let suitable partitions Q™ = {0 = s™ < s™ < • • • < s^„_;^ < = t] be given 

such that the diameters maxj=i A/m(s™ — ^^-i) 6° to zero. By exploiting again the stability 

(e(sjii), z(s;'!_i)) e ^(sjii) for j = 1, . . . , M™, we obtain that 

Wp(£(s;"), z(,s™)) - a(s;") : eisf) + i^(z(s™) - zisf_,)) 

> Wpieis]^,), zisf_,)) a(s™ J : £(s™ J - iaisf) - aisf_,)) : eisj) 

We shall take the sum above for j = 1, . . . , A'/™ and obtain that 

Wpieit), zit)) - ait) : £(t) + Dissi5(z, [0, t]) 

> Wpie^, zo) - aiO) : £0 - Y.^^'jisf) " '^(s" 1)) : eisf). (3.19) 
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Then, relation (|3.18p follows at once from Lebesgue dominated convergence since 



adr (s) : e(T'"(s))ds, 



where we used a standard notation for the piecewise mean on the partition Q™. In fact, eor™ and 
a dr converge to e and a at least almost everywhere, respectively, and e o r™ is uniformly 
bounded. Once (|3.18p is established, it is a standard matter to check that indeed Diss£)(2;, [0, t]) = 
4>{t) for all < e [0,T]. □ 

Finally, a consequence of p.l7p is the following Lipschitz regularity result. 

Corollary 3.3 (Lipschitz continuity). Under the assumptions of Theorem 13.21 if 
a G M^i'°°(0,T;M3^3), then we have (e,z) G M^i'°°(0, T; M^^^ x W^^J^). 

We shall complement the above detailed existence analysis by providing a local Lipschitz con- 
tinuous dependence result for the smooth case p > (see [43l Thm. 7.4]). 

Theorem 3.4 (Continuous dependence for p > 0). Let the assumptions of Theorem 13.21 hold 
p > 0, fTi, (72 G W^'^ {QjT;M.^^yj^) , suitably stable initial data (eo.ij-2o,i) '^'^'^ (^0,21^0,2) be given 
and (si,Zi) and (£2,2:2) be two corresponding energetic solutions to p.2p - (j3.3p . Then, there exists 
a positive constant c depending only on a, the bound and the Lipschitz constant of 'S/'^Wp, and 
||(Ti||^^i,i^Q j,.jj3x3-) for i = 1,2 such that 

|(ei-£2)(t)P + |(2l-22)(t)P 

< c (|eo,i " £0,2!' + 1^0,1 - zo,2f + ~ ^^2 ||^i.i(o,t;M-l,-^)) ^ P'^]- (3-20) 
Proof. Let us start by introducing some convenient notation. In particular, let 

ZiJ \OzWp[et,Zi)J 
V^M-,:=f^-^''('-"') a„I^p(£.2.A for, = 10. 



dezWp{ei,Zt) d,zWp{ei,Zi 

Next, by exploiting the above mentioned equivalence between p.2|) - (|3.3|l and ()3.4p . one readily 
checks that 

(VM^i - VM^2) ■ (yi - 2)2) <(ti-(T2): (£1-62) a.c. in(0,T), (3.21) 

where of course ■ is the scalar product in M^y,^ x M^^^. Moreover, we shall use £ := £1 — £2, 
z := zi — Z2 and so on. Within this proof, the symbol c will denote any positive constant possibly 
depending on a, ||W^pI|(72,i(k=*x-^ xr-^**-'')' ^^"^ ™ ll'''i|livi-i(o t-k^^'-^) * = 1,2. Let us define 

7 ■.= deW:e + dW:z > a\e\'^ + a\z\^ = a|y|^ 

where we also used the uniform convexity of Wp. Now, by differentiating 7 with respect to time 
and exploiting the smoothness of Wp, one gets that 

7 = (VVFi - VW2 + y^Wiy) ■ yi - {VWi - VW2 + V^WiV) ■ m 

< 2{WWi-VW2)-{yi-y2) 

+ I - S/Wi + VI^2 + y^Wiy\ \yi\ + \ - S/W2 + VM^i - y'^W2y\ M 

< 2a :^+c(|2)i| + |y2|)|f|' a.e. in(0,r). 
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By collecting the above computation we check that, for all t £ [0, T], 

7(0=7(0)+/ 7rfs<7(0) + 2a(t) :e(t)-2CT(0) :£o-2 / W:eds + cf i\y i\ + Ihl) 1 ds 
Jo Jo Jo 

< ^lit) + c [\eo\^ + \zo\^ + + Wm' + J\\yi\ + Ihlhds^ . 

The assertion follows by Gronwall's lemma. □ 

Properties of the approximations. The above detailed existence proof exploits a discrete 
construction which is interesting in itself. Let us condense in the following lemma the above 
proved results on the discrete scheme. Note that the result is less sharp for p — since we do not 
know whether the solutions are unique in this case. 

Lemma 3.5 (Convergence). Under the assumptions of Theorem 13.21 the incremental solutions 
(e",^") of problem (|3.7p for partitions P" with diameters r" going to zero are such that, possibly 
extracting a not relabeled subsequence, for all t G [0,r], 

z" — !■ z uniformly in [0,r], 
mssD{z",[0,t]) ^ DissD{z,[0,t]), 

£"(i) ^ sit), 

Wp(£"(t),z"(0)^W^p(£(t),z(0), 
for some pair {£,z) which solves p.2p - p.3p . As p>0 the whole sequence (£",z") converges. 

We conclude this section by recalling from [43] (see also [S^) an a priori error estimate of order 
1/2 for the above discussed discrete approximations. The latter error bound is however restricted 
the smooth situation p > 0. 

Lemma 3.6 (Error). Under the assumptions of Lemma 13.51 let p > 0. Then there exists a positive 
constant c depending on a, the bound and the Lipschitz constant of \7^Wp, || VFp||(^2,i(i{3x3 ^^3x3^, 

{eo,zo), and \W\\w^-'^{o,t-rI'}^) such that 

|(£-e")(t)| + |(z-z")(t)| <c(r")i/2 VtG[0,T]. (3.22) 

We shall not provide here a proof of the above lemma. Indeed, in case a £ W^^'°°(0, T; R^yj^) 
it suffices to rewrite in the current setting the argument of [37l Thm. 4.3]. Moreover, the proof 
can be adapted with little additional intricacy for the current absolutely continuous case a G 
PFi'i(0,r;R3^3) as well. 



4 Incremental minimization for the boundary 
value problem 

In this section we focus on a minimum problem which arises from the time incremental approxima- 
tion of the quasi-static evolution. Since we are actually dealing with a rate-independent evolution, 
this minimum problem is of course the basic tool for understanding the phenomenon. More- 
over, the study of the time discrete seems to be heavily addressed by the engineering community 
[24l[3Ql[3ll[32l|47l|48l[49j. Finally, the time incremental situation will turn out to be better suited 
than the time-continuous one in order to prove convergence of space approximations. 
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The data of the minimum problem are the current value z G L^{fl, I^dov) °^ ^^^^ inelastic strain 
and the updated values u^" e H^{^;M.^) of the boundary displacement and £ G {H^{n;M.^)Y of 
the total load. We shall be interested in solving the following 

{u,z)e ArgMin {Wp.^{v,w) - {£,v) +V{w -z)). (4.1) 

The existence of minimizers to the latter problem is a straightforward application of the Direct 
Method of the Calculus of Variations [TU] . Indeed, {v,w) i-^ Wp^^{v,w)+T>{w—z) — {£,v) is trivially 
coercive and lower semicontinuous with respect to the weak topology in and y'^{u^") is convex 
and closed. As far as uniqueness is concerned one should observe that Wp.^ is uniformly convex 
for all p, ly >Q. 

Let us state here a preliminary lemma whose proof can be obtained by means of standard 
computations on the quadratic form C. 

Lemma 4.1 (Change of boundary conditions). Let u^"\ v^" e H\n;M.^), z e L2(ri, R''^^^), and 
£ e (iJi(ri;M3))' he given. Moreover, let {u\z*) G J^''(u°") solve and v* = u* - vP" + v^'" . 
Then {v*,z*) solves 

(i;*,z*)e ArgMin (wpAv,z)+(c{e{v)^z):e{u^''^v^''^)~{£,v)+V{z-z)]. (4.2) 

On the other hand let (w*,z*) solve (g^l). Then {v* - v^" +u°'",z*) solves (gT]). 

Problem (|4.ip is Holder continuously stable with respect to perturbations on the data z, u^", 
and £. Indeed, we have the following generalization of Lemma l3. II 

Lemma 4.2 (Continuous dependence). Let p, i'>0 be fixed and zi, Z2 € L^(ri,R^^^), u^" , £ 
ffi(f7;E3), and £1, £2 G {H^{n;R^)Y be given. Moreover, let iu^,Zi) G yiuf'") solve gT]) with 
u^" = uf", J = Ji, and £ = £i for i = 1,2. Then, there exists a constant c depending on cq, a, 
Cd, and C such that 

\\U1 - '«2||h1(0;K3) + \\Z1 - 22||^2(0;R3^x3) + '^ll^l - ^2||^1(j^.]k3^x3-| 

< c (h?"- - u^"\\'m^nm + 11^1 -^2lLi(n;K-3) + \\£, ~ £2\\lHHn;M^)y) ■ (4-3) 

Proof. We simply adapt the argument of Lemma |3 . 1 1 Owing to the minimality of (wi, zi) and the 
uniform convexity of Wp,iy we readily deduce that, for any {vi,wi) G y'^{uY"), 

"IkK -«l)lli2(0;K^^x3) +«||Z1 -?i;i||^.(^.jj3^x^^ 

< )^p,Avi,Wi) - {£i,vi) + V{wi - zi) 
-Wp,^(wi, zi) + (£1, ui) - P(zi - zi). 

On the other hand, the minimality of (u2, Z2) entails that, for all {v2, W2) G y'^{u§"), 

< Wp,,{v2,W2) - {£2,V2)+V{W2-Z2) - Wp,,iu2, Z2) + {£2, U2) ~ V{z2 - Z2) . 

By choosing {vi,vui) ~ {112 — + iif", 22) and {v2, W2) = (wi — w?" + v^" ^ zi) and taking the 
sum of the corresponding inequalities one easily deduces that 

a\\e{ui - U2) ~ e{uY" - uY")\\l,^^,^s.^^+ a\^^^ 

< 2C(£(uf ' - w^"-)) - / C(£(ui - U2) - (zi - Z2)) : e{uY' - u^') 
+2P(zi - Z2) + {£1 - £2, ui - U2) - {£1 - £2,uY" - uY"). 
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Hence, we readily find a positive constant c depending on a, Cd, and C in such a way that 

||e(wi - '^2)||^2(n;R3X3) + ll^l - 22|li2(n;R3^x3) + '^ll^l - ^2||^1(0;R3^x3) 
< C - U°''||h1(M3) + 11^1 - ^2|lLl(a.R3^x3) + 114 - 4|l(ifl(M3))') • 

Whence, the assertion foUows from Korn's incquahty (|2.ip . □ 

Convergence of space approximations. Let us now turn our attention to some space ap- 
proximation procedure and recall the material of Section [21 We denote by 3^^q the set y^^Q := 
yj^ny'^{0). Given {u, z) — p'j[{u, z) wc shall also denote by pf^ i{u, z) := u and pf^ 2('^> ^) •= ^- For 
the sake of completeness, we shall consider also some approximate situation. Indeed, we ask that 
for each [u^" ,z) G and I G {H^{n-M?'))'), there exist {uf\zh) G 3^/^ and 4 G {H^{n-M.^))' 
such that 

(u^"-, z„) ^ 2) strongly in M^^) x L\n- W^l^), 

and 4-^^ strongly in {H^ {VL\W^))' . (4.4) 

We shall be concerned with the approximating minimum problem 

{uh,zh)e ArgMin {WpA^, z) - {Ih^u) ^ V{z -zh)) . (4.5) 

The latter problem is of course uniquely solvable since (u, z) i-^ W'p,i/(u, z) — {£h, u) + ^(z — z) is 
again uniformly convex, coercive, and lower semicontinuous in yf^ and yj^ q is convex and closed. 

Assuming (j4.4p and letting (u, z) and {ufi,zi^) solve the minimum problem (|4.ip and (|4.5p . 
respectively, the main issue of this section is that of proving that {uh,Zh) converges to {u,z) 
strongly in . More precisely, in the case p > 0, some quantitative error estimates can be 
obtained. 

Lemma 4.3 (Error for p > 0). Let p > 0, > be given and (u, z) and {uf^^ z^) solve (j4.ip 
and (|4.5p . respectively. Moreover, let 

(4,t'-Ph,i(i',w')) = for all {v,w)ey andh>0. (4.6) 
Then, there exists a positive constant c depending on p, cg, a, Cd, and C such that 

\\u - U/i|[Hi(f2;R3) + 11^ ^ ^''IIl2(J2:R3^^3-, + l^\\z — ^'i||^1(0;R3x3-) 
< C - Uh"\\H^n:M^) + \\'^'~Zh\\mn;RlxJ)) 

+c(lK-4||(ifi(0;R3)), + ||z-K,2(«:^)llLi(n;R3^^3)) . (4.7) 

Let us comment that (|4.6p turns out to be fulfilled in the frame of conforming finite elements. 
Considering for simplicity the case where p'j^ ^ does not depend on ui, a fairly usual choice for £h 
is 

(4,^;):=(4K,iW) Vz;GZi, 

whence (|4.6p follows. 

Proof. The estimate follows by carefully reconsidering the continuous dependence proof of Lemma 
14.21 and exploiting Galerkin's orthogonality (|2.6p . Indeed, making use of Lemma 14. 1[ one obtains 
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for V = u — u and Vh ^ — u)^" , 

- ^''«)lli2(jj.R3x3) + a\\z - Zh\\l,^^.^3.j^ + av\\z - Zh\\\^ ^^^.^...^ 

< A,{vh, Zh) + gp{zh) + / <C{e{vh) - zn) : e(u°"') + V{zh - z) - w^. - v) 



-A,{v, z) - Gpiz) - / Cisiv) - z) : - P(z - z) (4.8) 

where we have denoted by Gp : L'^{n,R^^^^) [0, +oo] the convex functional 



Gpiz) :=.Fp(z) -C2||z||^,(,j,j3^x3). 
Moreover, arguing exactly as in Lemma 14.21 and defining (-5, z) := p'j^iv, z), we readily check that 

< A(w, S) + Gpiz) + / C(£(^) - i) : e{uY") + V{z - zh) - (4, v ~ v,,) 



-A^.{vh, Zh) - Gp{zh) - / ^{e{vh) - zn) : £«"■) - V{zh - zh). (4.9) 

Taking the sum of the latter inequalities and exploiting ()2.7p . (14. 6p . and (it)^" , 0) G 3^/^, we easily 
check that 

a|k(« - ^'/«)lli2(a;R3x3) +a\\z- Zh\\l, ^^.^...^ + av\\z - Zh\\, 
< f C{e{vh -v)- {zh - z)) : e{u'''' - uf^ + 2V{z - Zh) 



+ (£ - 4, f - Vh) + Gpiz) - Gpiz) + Viz - S), 
and the assertion follows. □ 

We shall now turn to some (necessarily weaker) quantitative convergence estimate for the 
specific case p = 0. 

Lemma 4.4 (Convergence for p — 0). Under the assumptions of Lemma H.Sl let p = 0. Moreover, 
let iv,z) := Phiu—u^" , z) and iv,z) := (g/i(w— u^'''), r^(w)). Then, there exists a positive constant 
c depending on ci and the same constant of (|4.7p such that 

\\u - U/i|[hi(0;K3) + 11^ ~ ^''llL2(n:R3^^^) + '^W^ ^ ^'iIIhi(0;R=*j^^) 

< C (^\\U^" - ?^)?'1^1(0;E3x3) + IIZ - Z^||^i(f^.jj.3^x3) + \\e - 4||(^/l(n;R3))' + IN " ^ 1 1 L 1 (f2;R-7J ) ) 



+c i^A^iv, z) - A(w, S)+ J C(£(u - i)) - (z - S)) : eiuf^") 

+c((4,v-«) + ||5-z|Li(m3^x3)) . (4.10) 

Since of course phiv^w) — (g/i(w), rj^(u')) strongly converges to zero in , estimate (|4.10p 
proves in particular that, assuming (|4.4p . the strong convergence of the approximations holds. 

Proof. This proof follows the same lines of Lemma H31 We shall however replace (|4.9p as follows. 

< A(w, z) + cip||^i(f^.jj3^x3) + J Cieiv) - z) : £«'') + P(z - z,,) - {£,,,6 - vh) 

-A^ivh.Zh) - ci\\zh\\^i^^-i.;^3>c3,-^~ / Cieivh)-Zh) : eiu^") -Vizh-Zh), 



and again take its sum with (|4.8p . In order to reduce to the situation of Lemma [4.31 one needs to 
simply add and subtract the term z in most of the occurrences of z. This procedure of course 
produces the extra residual terms that appear in the last two lines of (|4.10|) . □ 
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5 The incremental problem. 

We shall prepare here some material in the direction of the full time-stepping procedure. To this 
aim, we assume to be given a partition P := {0 = to < ti < • • • < ijv-i < = T} with diameter 
T = max,=i,...,^(i, and data {uf"}l, e {H^{n;R^))^+\ {e,}^, e {{H\n;R^))T+\ 

and {uo,zq) G y'^{u^"). Hence, we find iteratively the unique solutions {{ui, Zi)}fLi to the 
problem 

{u„Zi)& ArgMin {Wp,^{u, z) - {£i,u) + V{z - z,^i)) for i ^ 1, . . . , N. (5.11) 

We shall denote by (m, z) the incremental solution which interpolates right-continuously the values 
(ui,Zi) on the partition P. Hence, the following a priori estimate holds true. 

Lemma 5.1 (A priori bounds). Let p, v > 0. Then there exists a positive constant c depending 
on a, yVp.jy{uo,zo), {£o,uo), and Y^iLi Ui ~ ^i-iWimin-M"))' that 

Wp,^(u, z) + Dissi,(2,[0,r]) < c. (5.12) 

Proof. From the minimality of (ui,Zi) in (jS.lip one has that 

>Vp,iy(Wi, Zi) - (li,Ui) + T){zi - Zj_i) 
< >Vp,t.(ui_i,Zi_i) - - {£i - ti-i,Ui-i). 

Taking the sum in the latter relation for i = 1, . . . , m, m < N , one has that 

m 

Wp,,.(u„i, Z„) - {Im, Um) + ^ 'D{zi - Zi_i) 

i=l 
m 

< Wp,i.(uo, 2o) - (4, Wo) - '^{^i - 

and the assertion follows from the uniform convexity of Wp.^ and the Gronwall lemma. □ 

Let us collect here some remark on the incremental problem (|5.11|1 in the space discretized 
situation. To this aim we shall refer to the notation introduced in Section[2]and assume to be given, 
for aU h>Q., suitable data {u^)l}to e {kh]to ^ iiHHn;R^))Y+\ and {uo,h,zo,h) 

such that (ito,ft, — u^^^, zg^h) G 3^o ft,- Hence, by solving iteratively the minimum problem, we define 
the right-continuous piecewise constant incremental solutions (uh, Zh). 

First of all, one should notice that the a priori bound of Lemma [5T] holds for (u/j, Zh) as well 
(of course the dependences of the constant arc referred to the approximating data) . Secondly, we 
are in the position of obtaining for {uh,Zh) the same continuous dependence as in Lemma 14.21 
This fact entails the convergence of the space approximated incremental problem in N steps to 
the corresponding limit. In particular, employing Lemma 14.31 or respectively, and performing 
an induction over i = 1, . . . , A^, we have the following result. 

Lemma 5.2 (Convergence for N steps as h — > 0). Under the above assumptions, let the pa- 
rameters >Q, N € N be fixed and assume that ufjl^ — > uf" in H^{n;M.'^), ii^h £i in 
{H^{n;R^)y, and {uqm^zqm) -> {uq,zo) in H\n;R^) x ^^(r^; M^^^^) as h ^ 0. Then, we have 
that Ui^h Ui in iJ^(f2;M^) as well, for all i = 1, . . . , N . 

Indeed, we would be in the position of stating a more precise quantitative bound for the error 
maxi<i<jv — Ui^;i||//i(0;R3) in terms of data. This bound will however deteriorate and eventually 
explode as A^ — > -|-cxd. 
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6 The evolution problem 

We shall finally turn to the study of the time-continuous problem. In particular, we are inter- 
ested in energetic solutions to (|1.5[) - (|1.9[) along with the above prescribed boundary displacement 
and boundary traction conditions. Namely, our solutions will be functions t i-^ {u{t),z{t)) G 
yiu^'^it)) such that t {£{t),u{t)) is integrable and, for all t G [0,T], 

iuit),z{t)) e {{u,z) G yiu^^it)) such that, V(?I,z) G y'iu^'^it)), 

Wp,Au,z)^{e{t),u) <WpAu,z) - {£{t),u) +V{z -z)}, (6.1) 
Wp,^(w(t), z{t)) - {£{t), u{t)) + Dissi,(z, [0, t]) 

= >Vp,,KO),z(0))-(^(0),w(0))- / {l{s),u{s))ds. (6.2) 

Jo 

Following the argument of Section [31 we arc in the position of proving the equivalence of the 
two formulations (|1.5p - (|1.9p and (|6.ip - (|6.2[) as soon as the above mentioned boundary condition 
(plus an extra homogeneous Neumann type condition for z when > 0) are considered and the 
solutions are assumed to be at least absolutely continuous. The latter is of course a quite natural 
regularity requirement and we will readily recover it in our framework. 

The main issue of this section is to fix i/ > and exploit the analysis of O |43] in order to 
obtain some existence, uniqueness, and convergence of approximations result. Apart from infinite 
dimensions, the arguments involved here are quite close to those of Section [3l Owing to this 
consideration, we will mainly sketch the proofs of the forthcoming results by heavily referring to 
the corresponding material in Section [31 



An equivalent problem. It is convenient to introduce yet another equivalent formulation of 
problem (|6.ip - (|6.2|) by replacing the variable u by v ~ u — u^^'-' . The main advantage of this 
change of variables is that the energetic formulation for (w, z) takes values in the fixed phase space 
:= y^(0). Indeed, in the same spirit of Lemma [4.1) one readily computes that 

WpA^, z) - {£, u) = WpAv, z)+ f C{e{v) - z) : e[u^'') - {£, v) + C[e{u^'')) - {£, v^"). 

Jn 

Hence, one checks that (m, z) is an energetic solution if and only if {v, z) : t t-^ y^ is such that, 
for all t G [0,T], 

{v{t),z{t)) G Sit) := |(w,2) G y^ such that, V(w,z) G J^q. 

WpAv, z) - {L{t), [v, z)) < WpAv,z) - {Lit), {v,z)) + Viz - z))}, (6.3) 
WpAvit),z{t)) ~~ {Lit), ivit), zit))) + q{t) + Dissi,(z, [0, t]) 
= WpAviO), z{0)) - (L(0), (viO), z(0))) + g(0) 

{l{s),vis))ds~ / {lis),u^"is))ds, (6.4) 



where we have denoted by L : [0, T] (3^0 )' the functional 

{L{t),iv,z)):=- [ Cieiv)~z):eiu''"it)) + {£it),v) yiv,z)eyi;, te[0,T]. 



Here (•, •) is used for the duality pairing between (3^o )' and 3^o i well. Moreover, the function 
(J : [0, T] -> R is defined as 

qit) := C(wD-(t)) - {eit),u°'^-{t)) yt G [0,T]. 
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We shall explicitly observe that u^'' G W^'^{0,T; H^{n;R^)) and £ e W^^^{0,T; {H^{n,R^)y) 
entail that L e W^'^{0,T; (y^;)') and qeW^'\0,T). 

From now on, we will focus on problem (|6.3p - (|6.4p and leave to the reader the straightforward 
interpretation of the forthcoming results for our original variable u. Let us start from the following 
existence result. 

Theorem 6.1 (Existence for v > 0). Let > and p > 0. Given L <E W^'\Q,T; {^{0))'), 
q G W^'^{0,T), and {vq, zq) G 5(0), there exists an energetic solution (w, z) to (|6.3p - (|6.4p such 
that iv{0),z{0)) = {vo,zo). Moreover {v,z) e W^'^{0,T;y^). 

We shall not provide here a full proof of this result. Indeed, it suffices to suitably adapt the 
machinery of Lemma lS^ to the situation of ()6.3|) - (j6.4|l . In particular, we argue again by discretizing 
the problem on a sequence of partitions P" with diameter going to zero. The corresponding 
incremental problems 

iv^,z,)eATgMm{Wp.,^{v,z)-{Litf),u)+Viz~z'^_,)) for i = l,...,N'\ (6.5) 

will turn out to be solvable by means of the results of Section S) Namely, we can introduce 
some right-continuous and piccewise constant intcrpolant (w",z") of the discrete solution on the 
partition P". Moreover, we exploit Lemma l5 . 1 1 which entails that 

sup Wp^u{v"{t), z"{t)) and Var[o 7-](z") are bounded independently of n. 

te[o,T] 

Indeed, the latter bound depends now on Wp,,y(uo, zg), ||-^||vFi.i(o.T;(i;''(o)))')' ^^'^ lkllwi i(o,T)- 

As for the limit, we will make use of some extended version of Holly's principle |28[ Thm. 3.1] 
and find a (not relabeled) subsequence of partitions and a non-decreasing function cj) : [0, T] — » 
[0, +oo) such that 

z"(t)-^z(t) weakly in H\n;R^^l^) and Dissp(z", [0, i]) ^ 0(i) for aU te[0,T], 
Dissi,(z, [s,t]) < (j)it) - (j){s) V[s,t] C [0,T]. 

Indeed, here we have used in a crucial way that v > 0, i.e., the sublevels of Wp^i, are compact 
in L'^{n;Rl^j^) x L^{n;R^^^^). Moreover, we have that v"(t) = Cz^'it), C being linear, and 
Cz'\t) Cz[t) = v{t) weakly in H^{n]R^) for all t G [0,r], where {v{t),{)) G y^. 

The set of stable trajectories S := Utt=[o,T]{tj<S{t)) is closed with respect to the weak topology 
of . Namely, letting [tk^Vk, z^) G S with tk ^ t and (vk,Zk) —> (w,z) weakly in 3^q, we 
readily exploit the lower semicontinuity of Wp^,y, the weak continuity of V in {fl;R'^J^^) , and 
the continuity of L and get that 

Wp.^iv, z) + {L{t), {v, z)) < liminf {WpAvk,Zk) + {L{tk), {vk,Zk))) 
< liminf (>Vp.,(w,z) + {L{tk),iv,z))+V{zk-z)) = Wp,Av,z) + {L{t),{v,z)) +V{z -z) 

k — ^+oo 

for all {v,z) G 3^Q . Namely, {t, v,z) £ S and the stabihty condition (|6.3[) easily follows. Moreover, 
the initial condition is fulfilled by construction and the uniform convexity of Wp,i/ along with 
stability entail that the whole sequence £(w"(i)) actually converges to e{v{t)). 

As for to prove that {v, z) fulfills (|6.4p we readily deduce from the above stated convergences 
and lower semicontinuity arguments (see p.l4p ) that the equivalent of p.lSp holds. Indeed we 
have that 

>Vp..(«"(0, z"(t)) - (i(r"(t)), (v^it), z^m + q{r-{t)) + Dissp(z", [0, r"(t)]) 
< Wp,,(wo, zo) - mO), («o, 2:0)) + q{0) 

{e{s),v"{s))ds- {e{s),u{s))ds. (6.6) 

Jo Jo 



21 




{e{s),u{s))ds. (6.7) 



Moreover, again by stability, one has that Wp,^ z^{t)) — > Wp.jy{v{t), z{t)) as weU (see p.l6p ). 
As a by-product, the above stated weak convergence for (t;"(f), z"(t)) turns out to be actually 
strong in y . 

Exactly as in Theorem 13. 2[ the absolute continuity of (w, z) follows at once from that of L 
and g, relation (j6.7p . the uniform convexity of Wp^i^, and stability (|6.3p . In particular, we are 
in the position of reproducing the same argument as in (|3.19p and, exploiting once more stability 
and the continuity of data, obtain the upper energy estimate as well. Namely, one has that 
(j){t) = Dissi5(z, [0,t]) for all t e [0,T]. The existence proof is hence complete. 

Again, energetic solutions corresponding to Lipschitz continuous data turn out to be Lipschitz 
continuous as well. 

Lemma 6.2 (Lipschitz continuity). Under the assumptions of Theorem \6.1\ whenever L e VF^'°°(0, T; (3^0 )') 
and q €W^'°^{0,T), we have {e, z) € W^'°°{0,T;y^). 

Existence by smoothness. The above sketched existence proof exploits in a crucial way the 
compactness of the sublevcls of Wp,,y for > in the weak topology of (Q; M.'^) x (Q; R ^cv) 
and works for any p > 0. An alternative approach to existence of solutions of the energetic 
formulation is however available in the smooth situation p > by means of the construction 
of [43l Sec. 7], for instance. A possible advantage of this perspective is that of gaining explicit 
convergence rates. We shall address this issue elsewhere. 

In the above mentioned smooth situation p > no compactness is assumed for energy-bounded 
states but the energy functional Wp^i, : y [0, +oo) is required to be C^'^. This again forces 
ly > 0. Namely, given h e C^'^{R) with h" e L°°(R), one has that the functional 
n : L^{n;M.f^^) M defined by 



is C^'^ if and only if h is quadratic (and in this case Ti G C°°). On the other hand, H is C^'^ 
on i?i(r2;R3^^^). This fact entails that Wp.^ is C^-i on y if and only if v > 0. 

Continuous dependence. We are in the position of reproducing the continuous dependence 
result of Section [3] in the present framework and for p, i> > 0. Once again continuous depen- 
dence relies on uniform convexity and C^'^ continuity of the energy functional. In particular, the 
assumption v > 0, which of course plays no role in Lemma [3^ is actually needed here (see above). 

Properties of the approximations. The time discretization technique described above has of 
course some interest in itself. Let us collect for convenience some related result in the following. 

Lemma 6.3. Let v > 0. Under the assumptions of Theorem \6.1\ the incremental solutions 
(?;",z") of problem (|6.5p for partitions P" with diameters r" going to are such that, possibly 
extracting a not relabeled subsequence, for all t G [0,T], 




z"^z strongly m C{[0,T]; H\n;Rf^^)), 
Dissi,(z", [0,t]) ^ Dissi,(z, [0,t]), 
v^it) -> v{t) strongly m H^{n;R^), 
WpAv"{t),z"{t))^WpAv{t),z{t)), 
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for some (v, z) which solves (|6.3p - (|6.4p . As p > the whole sequence is convergent to the unique 
energetic solution {v,z) and there exists a positive constant c depending on a, \\'Wp,v\\c^.'^(y<^-;K), 
(wo,zo), \\L\\w^-^(o,T:iy(a))'), and ||q|| wi.i(o,T) such that 

\\{v ^'")(^)||HMO;E3) + \\{z ^")WllH^(n;K-3) < c{T-f'' Vt G [0,T]. (6.8) 

Full space-time approximations. We conclude this analysis by commenting on the possibility 
of performing a full space-time approximation of the problem. To this aim let us refer to the above 
introduced notations, consider some approximation parameter h > 0, and reduce the energetic 
formulation ()6.3p - (|6.4p to the spaces y^^Q exhausting J^q- We shall be considering in particular 
some discrete values i, i)}ilo defined inductively from suitable initial data {vh,o,Zhfi) & 
letting (uq , Zq) = (uh.Oi -s^/i.o) a-nd solving the following incremental problem 

e ArgMin(>Vp,.(«,z)-(L(^n,«)+2?(^-<.-i)) for z = l,...,N^. (6.9) 

Again, the unique solvability of the latter problems is ensured by uniform convexity and lower 
semicontinuity, i.e., it is independent of h. We will denote as usual by [v^, z^) the corresponding 
incremental solutions. 

Our first observation is that, arguing exactly as above, whenever the assumptions of Theorem 
16.11 are fulfilled and the initial data are bounded in energy independently of ft,, the usual bound 

sup W p^u{vh{i) Zh{t)) and Diss-p(2;JJ, [0, T]) arc bounded indep. of n and h, (6.10) 
«e[o,T] 

can be obtained. 

Convergence for the space-discretized problem. Assume ft, > 0. Then, we are in the 
position of reproducing the argument of Theorem 16 . 1 1 and deduce the existence of a limiting space- 
approximated energetic solution [vy^ Zh)- To this aim, the restriction v > Q could even be avoided 
whenever are chosen to be finite dimensional, for instance. Moreover, the fully discrete solution 
{v^, z^) converges to (w^, Zh) in the sense of Lemma 1673) as n +oo. We shall not give a detailed 
proof of these facts but rather limit ourselves in observing that the energetic formulation (|6.3|l - (j6.4p 
can be rewritten in q with no intricacy. In particular, estimate ()6.10p is again the starting point 
for the limit procedure. 

Once the energetic solution {vh,Zh) ■ [0,T] yji^Q is found (uniqueness again follows in case 
p > 0) we are in the condition of considering the limit as ft goes to zero as well. To this aim, we 
shall assume that the corresponding initial data converge together with their energies, namely 

Wp,,ivh,o,Zh,o) - (i(0), {vh{0), zh.(0))) ^ Wp,,(wo, zo) - (L(0), («o, zq)). 

In this case, it is straightforward to check that the bound ()6.10|) is preserved while passing to the 
limit in ft. Assuming > 0, this entails the possibility of extracting a (not relabeled) subsequence 
pointwise converging to an energetic solution {v, z) : [0, T] . In case p > 0, the latter is indeed 

the unique energetic solution whose existence is stated in Theorem 16. II In order to check this we 
briefiy comment on relations (|6.3[) - (|6.4p . As for (|6.3p . let us fix t £ [0, T] and any (w, z) G 3^o ^^'^ 
exploit the stability of {vii{t), Zh{t)) in order to get that, for all (iJ, z) G , 

WpAvhit),Zh{t)) - {L{t),{vhit),Zh{t))) 
< WpApI{v,z)) - {L{t),pl{v,z)))+'D{^^h-plA^,z)). 

Hence, the stability of {v{t),z{t)) follows by passing to the limit in ft. As for the upper energy 
estimate we fix a uniform partition := {s™, j = 0, . . . ,M : s™ = jt/m}, exploit the upper 
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energy estimate for (vh,Zh), and get that 

m 

WpAMt),^^h{t)) - {L{t),{vu{t),zu{t))) + q{t) +Y,V{zu{sf) ~ zu{sf_,)) 

< y^pA^hfl.Zhfi) - (i(0), {vhfi, Zhfl)) + g(0) 

{t{s),VH{s))ds- [ {i{s)Ms))ds. 
Jo 

It hence suffices to pass to the hmit in h first and then in m in order to get the upper energy 
estimate for {v,z). Finally, the lower energy estimate for {v,z) follows as above from the upper 
energy estimate, stability, uniform convexity of yVp,^, and the continuity of L and q. We refer 
to [ini m] for a full proof of the above convergence argument. However, we shall remark that no 
quantitative estimates for the approximations arc given. 

Convergence for the time-discretized problem. Let us consider now the limit as h goes 
to first. Owing to Lemma 15.21 we are in the position of establishing a (quantitative) strong 
convergence result for the corresponding time discretized solutions (w",z"). Indeed, one could 
exhibit some explicit error control which however explodes with n. Moreover, in the case v > 0, 
since (w", z") are uniquely determined, the subsequent limit in n can be taken exactly as above 
and the convergence to an energetic solution (v, z) is ensured. 

Joint convergence. Assume now u > 0. Owing to (|6.10[) we are of course in the position of 
passing to the limit with respect to both n and h simultaneously in (wJJ, zjj;). By arguing as above 
the stability of the limit (w, z) will follow at once by using the closedness of S and the convergence 
of projections. As for the upper energy estimate, we combine the above exploited techniques and 
pass to the liminf in the following relation (see (|6.6p ) 

>Vp,.(^;^(^),4(^))-(L(r"(^)),(^;;:(^),z;,'(<)))+g(r"(^)) + Dissp(z;^,[o,r"(^)]) 

< yVpAvo,h,zo,h) - (^(0), ivo,h, zoA) + 9(0) 

- r ^'\e{.s),v]i{s))ds- r ^'\ii.s),ui-'^))ds. (6.11) 
Jo Jo 

Once the upper energy estimate is established, the uniform convexity of Wp.i, the continuity of L 
and q, and the stability of {v, z) entail that also the lower energy estimate holds. Namely, (w, z) 
is an energetic solution to (|6.3p - (|6.4p and it is unique as p > 0. 

Of course, whenever p > we would be able to show some convergence of order 1/2 in time. 
On the other hand, by passing to the limit in time we loose the chance to estimate the error in 
space (see above). Hence, so far we are not able to provide an explicit space-time error bound for 
the joint limit procedure. 

7 The limits p, ^ 0. 

Up to this point, the parameters p and u have been systematically assumed to be fixed throughout 
the analysis. The limit — > is however of some interest since it describes the behavior of the 
model toward its non- regularized limit. As for p wc have to mention that our modeling choice 
corresponds to the limit situation p ~ Q . On the other hand the smooth situation p > is better 
suited for numerical implementation. Moreover, all problems arc continuously dependent on data 
for p > while energetic evolutions are not known to be unique for p = 0. 
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In this section wc shall discuss the possibility of obtaining suitable asymptotic results for p and 
(possibly) v going to zero within the constitutive relation, the minimum problem, the incremental 
problem, and the evolution problem. We will explicitly treat the space approximated case and 
discuss joint limits of parameters and time and/or space approximations. 

As a general remark, one should notice that the choice p = v ^ Q does not affect the wcll- 
posedness of the minimum problems since the uniform convexity of the corresponding functionals is 
preserved, this being true also for space approximations. Secondly, a priori bounds on sequences of 
solutions (either minimizing, incremental, or energetic) are usually available independently of the 
parameters. Whenever the compactness of sequences of solutions is obtained, the crucial feature 
in order to identify the limit of some possibly extracted subsequence is the F-convergence (see 
below) of the approximating functionals Wp (in the zero-dimensional case) and Wp,i/ (in three 
dimensions). 

F-convergence issues. Let us collect here some preliminary remarks on the convergence prop- 
erties of functions and functionals under consideration. The basic notion in this direction is of 
course that of F-convergence [18l [T9|. The reader is referred to the monographs [3l [11] for a 
comprehensive discussion. Let us however recall here that, given a metric space X and functions 
Qm g : X (—00, +00], we say that —>■ g in the sense of F— convergence in X iff 

g{x) < liminf (7„(a;„) Va;„ — > x and (''•12) 

n — ^ + 00 

Vx e X there exists a;„ x such that g{x) > limsup (7„(a;„). (7-13) 

We shall classically refer to (|7.12|) as T -liminf inequality and to x„ in (|7.13p as the recovery sequence 
for X. Moreover, letting X be a Banach space, we say that gn g in the sense of Mosco [3] if 
gn ^ g in the sense of F-convergence with respect to both the strong and weak topology of X. 

Let us mention that the issue of the convergence of rate-independent evolution problems un- 
der approximation is indeed a crucial one. A general abstract theory of F-convergence for rate- 
independent systems is detailed in |41| . 

Henceforth, we shall refer to the current choice (|2.3p and explicitly ask the function / to be 
convex and non-decreasing. This entails in particular that Fp F pointwise and non-decreasing. 
The smoothness of Fp and the latter convergence entail by means of O Thm. 2.40, p. 198] that 
Fp ^ F in the sense of F-convergence in K^^.^. As a consequence and by using [31 Thm. 2.15, p. 
138], we have that 

Wp -> Wo in the sense of F-convergence in M^y^j x M^^^^. (7.14) 

As for the three-dimensional situation, let us start by observing that J-p J- in the sense of 
F-convergence with respect to both the strong and the weak topology in L^(r2;R'^^^) (namely, 
J-p converges to in the sense of Mosco [3]). This fact follows at once from [3l Thm. 2.40, 
p. 198] and the convexity of J^p. For all > fixed, we readily deduce in a quite similar way 
that J-^p^i, converges to ^0.1/ in the sense of Mosco in H^{n;R^^^^). Let us make precise the latter 
statement with the following. 

Lemma 7.1 (F-convergence of the inelastic energy). Let Pk P > and — > j/ > be 

non-increasing. Then J-p^_uk ~^^p,v in the sense of Mosco in H^^'^^{Q]W^^^). 

Proof. The above discussion may be readily extended in order to cover the case Vk v > 0. Let 
us turn to the situation i> = and Vk > instead. Of course, the F— liminf inequality ()7.12p 
easily follows from the F-convergcncc J-p^, — > J-p and lower semicontinuity considerations. As for 
the recovery sequence, letting z G L^(f2;M^^^) be fixed, we shall define Zk as the unique solution 
to the singular perturbation problem 

Zk + iykJzk = z in {H\n;RfJ)y, 
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where J : I^d^v) ^ (^H^^; I^dcv))' the Ricsz map. Wc have that (see, e.g., Lions [57]) 

Zk^z strongly in 1^3x3) ^/"iVz^p^O. 

Moreover, whenever \z\ < C3 almost everywhere in fl, the same bound holds for all Zk by the 
maximum principle. Hence, we readily check that 

and the assertion follows. □ 

We shall now turn our attention to the convergence of stored energies and state the following. 

Lemma 7.2 (F-convergence of the stored energy). Let Pk P > and ^ v > Q he non- 
increasing. Then 'Wp^^Vk ~^ ^p,^ the sense of Mosco in . 

We will not provide the reader with a detailed proof. Of course, the argument can be easily 
reproduced by arguing along the lines of the proof of Lemma [77T] 

7.1 Constitutive relation 

Let us denote by (e, z)p,r the incremental solution to the constitutive relation on the partition 
P := {0 = to < ti < ■ ■ ■ < In-i < tN ~ T} with diameter t, namely the right-continuous 
piecewise constant interpolant on the time partition of the solutions {(Sp, z^)} to 

(e;,z;)e ArgMin {Wp{e, z) - ait,) : e + Diz - z^-')) i = l,...,N, 

where a £ W^^'^(0, T; Rgy,^) and {s°,Zp) = (eo,^o) are given. Moreover, for aU p > 0, we wiU 
denote by (e, 2)p,o a solution for the time-continuous constitutive relation. Of course we would be 
in the position of considering approximating data Cp and (eq, ^o)p,t as well. We limit ourselves 
to the above situation just for the sake of simplicity. The main result of this subsection is the 
following. 

Theorem 7.3 (Convergence for the constitutive relation). Let Pk ^ P > and ^ r > 

either being constant or converging to 0. Then, possibly up to the extraction of a subsequence in 
the case {p,t) = (0,0), we have that 

(e, z)(p t-'ij, ^ (e, z)p,T pointwise in [0, T]. 

Indeed much more is true since the convergence of the component ^(p is uniform and we have 
convergences also of energies and dissipations. Moreover, one could consider the limits pk p > 
and/or rj, — > t > as well (which we however believe to be less interesting). We limit ourselves 
to the above statement for the sake of clarity. 

The situation of Theorem 17.31 is described in Figure 1 below where every parameter choice 
(p, r) in the p x t square gives rise to a solution either of the incremental problem (for r > 0) or 
the time-continuous problem (r = 0). Of course this solution is unique if (p, r) 7^ (0,0). Theorem 
17.31 entails that all the depicted limits (arrows) can be performed. 

Proof. By referring to Figure [1] we shall proceed by discussing limits of type a, b, c, and d. 

Limits of type a, namely {p,T)k —>■ (p, 0). These limits follow directly from Theorem 13.21 

Limits of type b, namely {p,T)k ^ (0,t) with r > 0. Since the time partition is fixed, the 
convergence of the whole sequence (e, z)(pj,^^') to the corresponding incremental solution (e, z)(o_t-) 
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Figure 1: Convergences for the constitutive relation 



is ensured by the r-convcrgcncc of the corresponding energy functionals. their equi-coercivity with 
respect to p, the continuity of R, and the continuous dependence of the incremental problem for 
p>0. 

The limit c, namely (p, r)fe s- (0, 0). Let us now turn to the joint limit. Again, the usual energy 
and dissipation bounds may be obtained and, by suitably choosing not relabeled subsequences, we 



find (£,z) : [0,r] ^ xR3^^3 g^chthat ^ z{t) and e(p,.),(t) ^ e{t) for all t e [0,T]. 



As for to prove the stability of {e{t),z{t)) we simply need to specialize the closure argument in 
Theorem 13.21 bv considering the parameter dependence on p. Here, the F— convergence (|7.14p is 
again crucial. In particular, let us redefine (see ()3.ip ). for all p > 0, 



and Sp := Ut£io,T]{t, Sp{t)). Owing to the F— convergence (|7.14p and the continuity of a we 
readily check that, for all {tp,£p,Zp) ^ Sp such that (tp,ep,Zp) converges to (to,£o>^o) as p 
one has that (tQ^eo,zo) G 5*0. As for the upper energy estimate, we readily pass to the liminf in 
the discrete upper equality estimate (|3.14p by means of the F-convergence (|7.14p and the fact that 
Wp — *■ Wq pointwise. Finally, the full energy equality follows again from stability. 

The limit d, namely {p,0)k — * (0,0). We shall not discuss this limit in detail since it follows 
easily along the lines of limit c above. □ 

7.2 The minimum problem 

We investigate for simplicity the situation of fixed data e H^{n;R^), g e (iJ^ (rj; M^))'^ 

and z e i^(f2; Kjov)- course, some more general situation of parameter-dependent data could 
be considered as well (see also the forthcoming Lemma f7.6p . Moreover, let us introduce for the 
purposes of this section the notation Tp^^ : — > (— oo, +oo] as 



for all p, > 0. Problem (j4.ip has a unique solution {u,z)p^^ £ y^iu^") for all given parameters 
p, J/ > 0. Moreover, we readily check that yVp^y((w, z)p^^) turns out to be bounded independently 
of p and v. Hence, {u,z)p^i, is weakly precompact in . 

Moreover, we shall consider the space approximated situation described by the mesh-size h > Q. 
For the sake of notational simplicity, we reduce ourselves to the oversimplified situation of data 
independent of /i . In particular, we assume uP" £ Uh for h small enough and define yj^yP") ■= 



Spit) {(£, z) e M3^3 X Rf^l such that, V{e,z) e K^^^ x W 
Wp{e,z)^a{t) : e < Wp(s,z)~a{t) : e + £"(2 - z)}, 



3x3 
dov' 



(7.15) 



Ip,Au, z) := Wp^^iu, z) - {£, u) + V{z - z) V(m, z) e y 
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g + {u^" , 0) . As for the general case, the foUowing discussion has to be restricted to the situation 
where convergence (|4.4|) holds for the approximating data u^", £h, and z/j. Consequently, we will 
make use of the notation 

Ip.ujiiu, z) := Tp^^{u, z) for [u^z)^y][ and +00 otherwise in . 

We shall start by providing the following convergence result. 

Lemma 7.4 (F-convergence of Tp,u,h)- Let pk ^ p > 0, Vk ^ > 0, and h > 0. Then 

Tpi^Mk ^Ip.u in the sense of Mosco in , (7-16) 
1pk.i^k,h ^ ^p,u,h in the sense of Mosco in y^. (7-17) 

Moreover, let /ifc — > 0. Then 

'^(p,u.h)k ~* ^p.i' in the sense of Mosco in y^ . (7-18) 

Proof. The convergence in (|7.16p follows directly from Lemma [721 and the strong continuity of T) 

in L\n,m.^^^l). 

Convergence (|7.17p is also straightforward. Namely, the liminf inequality for weakly converging 
sequences is immediate and the construction of recovery sequences follows at once from pointwise 
convergence (recall that y^^ = hence no singular perturbation is needed here). 

The full convergence situation of (|7.18p deserves some comment. Given any {u,z) G y^ , we 
define 

{u,z)i^p^^^h)k ■= ('J/^kH^CJ^))- 
Owing to the convergence and boundedness properties of the projectors qh^, and r)^^ (see Section 
[1]), we readily deduce that (u, z)(^p ,^j^)i^ (w, z) strongly in and 

W(p,i.,/i)fc((^i,^)(p,i.,/oJ y^p,,^{'n,z). 
The lim inf inequality follows once again from lower semicontinuity. □ 

The main result of this subsection concerns the possibility of considering (possibly joint) limits 
in the parameters p, v, and h and is graphically represented in Figure 2 below. 



V 




(0,0,0) p 

Figure 2: Convergences for the minimum problem 

Theorem 7.5 (Convergence for the minimum problem). Let Pk ^ P ^ ^, i^k ^ ^ ^, and 

hk ^ h > either being constant or converging to 0. Then 

{u,z)(p^^j^)^ ^ {u,z)p^y,h weakly in y" (X if h > 0). 

This result, whose proof is not reported, follows at once from Lemma [7.4l and the equi-coercivity 
and uniform convexity of the functionals. The limits {p, v, h) — > (p, v, 0) where already discussed 
in detail in Section [D 
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7.3 The incremental problem 

We shall extend the latter asymptotics for the minimum problem to the situation of the incremental 
problem on the fixed partition P := {0 = to < ti < ■ ■ ■ < t^^i < = T}. To this aim let the 
data {u^"''}'iLo, {^'}iIo and the initial datum {u°,z°) be suitably given independently of p and 
V (for simplicity). Then, for all p, j/ > we are entitled to solve the incremental problem and find 
a solution vector {(Up.,y, 2p.y)}i^o- Now, arguing as above, we easily obtain that yVp,u{u^p.„, Zp.„) is 
bounded independently of p, v, and i. For all given p, > 0, i = 1, . . . , TV, and z G L^{^] I^dcv)' 
we introduce the functionals Jp,y{-, ■,'z) : y (— oo, +(»] as 

j;^,{u, z,z) ~ Wp,,(u, z) - {i\u) + V{z - z) V(w, z) e 

Moreover, possibly taking into account the space-approximated situation, one would need to in- 
troduce space approximated data {wlj^'^'^liLo' {^Lli^o ^'^^ the initial datum Let us 
however restrict ourselves to the (over)simplified situation where the latter can be assumed to be 
independent of h. For all p, i/ > 0, ft- > 0, i = 1, . . . , A^, and z g L^{Q.; K^g^), we shall make use 
of the functionals Jp^^^h{-^ -,2) : y {—00, +00] defined as 

Jp jj j^{u,z,z):=J'p^^{u,z,z) if (ti, z)e3^^ and +00 otherwise. 

Let us start from the following F-convcrgcncc result. 

Lemma 7.6 (F-convergcnce of J^p^i,^h)- Let pk ^ p > 0, i^k > 0, and h > 0. Moreover, let 
Zk^z strongly in 2.1(0; M^^^^). Then, for all i = l,...,N, 

Jp^^i,i^{-, ■,'zk) ^ J'p.i,{-, ■,'z) in the sense of Mosco in y", (7-19) 

^p\,^.j>(-: -.^fe) ^ Jp.r.,hi-, the sense of Mosco m yj^. (7.20) 

Moreover, let /i^ — > 0. Then, for all i = 1, . . . , N , 

^(P,y,h)^{-^ -^Zk) JlA-^ •,^) the sense of Mosco in y" . (7.21) 

We are not reporting here the proof of the latter lemma for the sake of brevity. Indeed, the 
argument may be easily adapted from that of Lemma l7.4l bv exploiting the strong continuity of T> 
in Li(r2;R^g^), its lower semicontinuity in L^(ri;M^^^), and the triangle inequality 1^^ . 

By using Lemma l7.6l and denoting by (ti, z)p,i, and (u, z)p^u.h the incremental solutions related 
to the parameter choice (p, v) and, possibly, the space approximation, the main result of this 
subsection reads as follows. 

Theorem 7.7 (Convergence for the incremental problem for u > 0). Let v > he fixed and 
Pk — * P, o-nd hk —>■ h > either being constant of converging to 0. Then, for all t G [0, T], 

{u{t),z{t))p^^i,,h^ ^ {u{t),z{t))p^y,h strongly in y" . 

Of course, we would be in the position of considering the case ^'fc — > i', Pk ^ P > 0, and/or 
hk h > as well. We however restrict to the above situation for the sake of clarity. 

Lemma 17.61 entails the convergence of the incremental solutions as soon as the strong con- 
vergence of Zp^^i, or Zp^^^^hk in -^^(^^; K'^ev) ensured. In order to obtain the latter from the 
boundedness of energy through compactness we are forced once again to restrict our attention to 
the case v > 0. The proof of Theorem 17.71 follows then by simply taking steps in i. 
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7.4 The evolution problem 

Owing to the latter discussion on the incremental problem (see Lemma I7.6p . we shall restrict 
ourselves to the situation > from the very beginning (note that existence is not known for 
V = 0). For all p, /i > 0, let us denote by [v,z)p : [0,T] ^ and [v,z)p^h ■ [0,r] y^j^ the 
solutions to the corresponding energetic formulations for h = and h > (here and in what 
follows we have assumed the data L, q, and the initial datum (u°, z°) to be fixed independently of 
all approximations). The latter solutions are known to exists and turn out to be unique for p > 0. 
Moreover, let (v, z)p^r and (u, z)p^T-,h denote the unique incremental solutions to the problem on 
a given partition with diameter r. 

A variety of convergence results for (w,z)p, {v,z)p^h, {v,z)p^r, and {v,z)p^r,h have already 
been obtained. This subsection will complement the above discussions and complete the picture 
of convergence results for the time-continuous evolution problem. In particular, as soon as > 
is fixed, we are entitled to take (possibly joint) limits in (p, r, h) as it is graphically depicted in 
Figure 3 below. 
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Figure 3: Convergences for the evolution problem {y > 0) 



The main result of this subsection reads as follows. 

Theorem 7.8 (Convergence for the evolution problem for v > 0). Let ly > Q be fixed and pk p, 
Tk ^ T > 0, and hk —f h > Q either being constant of converging to 0. Then, possibly extracting 
not-relabeled subsequences if {p,t) = (0,0), for all t G [0, T], 

{v{t),z{t))(^p^^^h)k ^ {v{t),z{t))p,r.h strongly in y^ . 

Sketch of the proof. Referring to Figure 3, let us start by observing that the limits of type a and 
b were already obtained in Theorem 17.71 and Theorem 16. 1[ respectively. Moreover, the limits of 
type c have been discussed at the end of Section [6l 

Limits of type d. This limits can be established by simply adapting to the current three- 
dimensional situation the argument of Theorem 17.31 In case h > 0, the latter adaptation is even 
simplified by finite-dimensionality and the convergence result would hold for v = as well. 

The limit e. By suitably extracting (not-relabeled) subsequences we readily find {v, z) : [0, T] — > 
y^ such that, for all i S [0, T], 

{v{t),z{t))(^p,r.h)k ^ {v{t),z{t)) weakly in y^ , 

z^p,r.h),{t) zit) strongly in L2(f};R3^^^3)_ (7 22) 

Hence, we are left to prove that indeed (w, z) is a solution of the evolution problem, i.e., check for 
the stability condition (|6.ip and the energy equality (|6.2p . 
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As for the former, we exploit Lemma 17.21 and, for all {v,z) G y^, by letting {v,z)k '■= 
{qh^ iv), rf^^ (z)) we check that 

WoAv{t).z{t))~{L{t).{v{i).z{t))) 

< liminf (Wp,^^{(v{tr^),z{tr,))(p,r,h)J ^ {HtrJ,{{v{trJ,z{trJ)(p,r.h)j) 

A:— >+oc \ / 

- ^i^l^^ (^Pk:i'{{v,z)k) - (i(irj, {{V,z)k) + V{zh - ^;(p,r,h.) j) 

/C — ^~pOO \ / 

= Wo,,(F, z) - (L(t), {v,z)) + V{z - z{t)) 

where we used some obvious notation for the point on the time-partition of diameter Tk such 
that < t — < Tfc, Lemma 17.21 the stability of (w, z)(p,T,/i)fc a-t time trt, and the strong 
continuity of V in L"^ [VL^W^^^I) . 

The upper energy estimate (and hence (|6.2[l ) follows by simply passing to the liminf as 
(p, T, h)k (0, 0, 0) in the discrete upper energy estimate (|6.6p . 

The limit f. This limit can be obtained along the same lines of limit e above, the argument 
being even simplified by the fact that here = and the upper energy estimate follows by passing 
to the liminf as (p,h)k — ^ (0,0) in the time-continuous upper energy estimate (|6.7p . □ 
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